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Abstract 



shock velocity when deuterium transforms efficiently into HD molecules which then dominate 
gas cooling. Above this velocity the shocked gas is able to cool down to the temperature 



1 Introduction 



We study the effects of HD molecules on thermochemical evolution of the primordial gas 

> ■ 

' behind shock waves, possibly arised in the process of galaxy formation. Wc find the critical 

o ■ 
\o 
l> 
o 

■ of the cosmic microwave background. Under these conditions the corresponding Jeans mass 

Qh| molecules heat shocked gas, and at larger redshift their contribution to thermal evolution 

Q , becomes negligible. 

> 

■ Formation of stars is intimately connected with the ability of gas to cool. In the metal-free 
primordial medium the radiative cooling is mainly provided by molecular hydrogen and its 
isotope analogue HD. In the expanding universe H2 and HD form after the epoch of recom- 
bination [1-5]. Due to a non-zero dipole moment and lower excitation energy HD molecules can 
efficiently cool gas at temeperature T < 200 K, where the rate of H2 cooling decreases sharply. 
If HD abundance is small, cooling of primordial gas stops practically at T < 200 K. 

Thus, thermal evolution of the primordial gas at low temperature, and as a consequence char- 
acteristics of the first stars in the universe critically depend on the abundance of HD molecules. 
At present, an analysis of conditions when HD can be thermodynamically important in the pri- 
mordial gas is absent. Conclusions about the role of HD molecules in a prestellar universe are 
contradictory. In particular, it is shown in jj] that HD cooling never dominates in a collapsing 
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spherical cloud. At the same time the authors [6-8] point out that HD cooling can be important 
in primordial clouds, although their calculations are restricted only by initial stages of the col- 
lapse. It is obvious though, that the abundance of HD crucially depends on thermal evolution 
of gas in the temperature range T > 500 K. The apparent contradiction about the role of HD is 
connected with differences in the initial conditions used in 1_ from one side and in [6-8] from the 
other. In addition, the H2 cooling function adopted in 1 is overestimated. Therefore, in order 
to make firm conclusions about the role of HD molecules in a prestellar universe it is necessary 
to investigate the efficiency of HD formation in a wider range of initial conditions. 

It is known that formation of H2 and HD molecules is largely boosted behind shock waves 
[10-13]. It is connected mostly with the fact that temperature and fractional ionization behind 
the shock waves increase, and as a consequence the rates of molecular reactions are enhanced. In 
a postshock gas cooled down to temperature ~ 10^ the fraction of electrons remains sufficiently 
high, X }^ 0.001, which favours rapid formation of H2 molecules in the catalystic reactions with 
H~ ions. At such conditions H2 fraction can reach ~ 10~^. Further cooling is mainly provided 
by H2 molecules efficient in the temperature range 200 — 7000 K. When lower temperatures are 
reached, T < 200 K, deuterium begins to convert into HD molecules due to chemical fractionation 
|14|ll5j. The contribution of HD cooling in energy losses increases when temperature decreases, 
and if it becomes dominant the gas temperature can fall down to several tens of degrees. One 
can expect that at least in a restricted range of shock parameters formation of HD molecules is 
as efficient as to provide such a predominance. We aim to study this possibility. 

In Section 2 we describe a thermochemical model of a shocked gas; the results are presented in 
Section 3; in Section 4 discussion and in Section 5 summary are given. Throughout the paper we 
assume a ACDM cosmology with the parameters (Oq, J^a, ilm, ^b, h) = (1.0, 0.71, 0.29, 0.047, 0.72) 
and deuterium abundance 2.62 x 10~^ |16j . 

2 Molecular kinetics behind shock waves 

In the absence of thermal conductivity and diffusion the thermochemical evolution of gas behind 
a shock wave can be described by a system of ordinary differential equations for a Lagrangian 
fluid element, which include equations of chemical kinetics 

Xi = F{xi,T,n) - D{xi,T,n), (1) 

and energy equation 

T = lY.i^iix^,T,n) - A,{x„T,n)] + l-n, (2) 

2 



where Xi is the fraction of i-th species, Fi{x,T,n), Di{x,T,n) are the corresponding formation 
and destruction rates, rj(x, T, n), Aj(x, T, n) are the heating and coohng rates. Chemical kinetics 
of primordial gas include the following main species: H, H"*", H~, He, He^, He^"'', H2, H^, D, D"*", 
D^, HD, HD+, e. The rates for collisional and radiative processes are taken from and for 

ion from 5 . The energy equation accounts cooling processes connected with H, He, He"*", 
He'^^, such as collisional excitation and ionization, bremsstrahlung radiation, recombination, 
dielectronic recombination, molecular cooling by H2 and HD, and Compton cooling. In the 
absence of external ionizing radiation the abundances of chemical species and radiative cooling 
are determined by collisions. Thus, the right-hand side of (1) and first term of the right-hand 
side of (2) are proportional to the gas density, and for convenience we can introduce the fluence 
r] through 



where t is time, n is number density; below the results are presented as functions of the fluence. 

For the H2 cooling functions we adopted the expressions given in ^7] , the HD cooling is taken 
from ^S], the other rates are from ^5]. In addition, the effects of interaction between molecules 
and the CMB photons are accounted [21 1151 l2Uj which imply that when gas temperature is close 
to the CMB temperature, H2 and HD molecules are populated by the CMB photons, and then 
collisionally transfer the energy to kinetic temperature, thus heating the gas. Therefore, gas 
cannot cool lower than the CMB temperature. We assume that the gas temperature behind the 
front jumps and reach the value 



where = 3/16 is for a shock propagating in a static gas, and = 1/3 is for a shock wave 
formed in a head-on collision of two flows (clouds) with equal velocities Vc j21j . Neglecting 
thermal conduction the evolution of each Lagrangian volume of gas behind the shock is isobaric 
[ini [T^ H2] I so that the density is described by 



where ^ is the average molecular weight. 



dri = ndt 



(3) 




(4) 



3 Formation of HD behind shock waves 



In the contemporary scenarious of structure formation in the universe the first protogalaxies 
form at the epoch z = 10 — 30. For specifity, we consider thermochemical evolution of baryons 
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behind a shock at the redshift z = 20. Formation of dark haloes (future protogalaxies) and their 
subsequent viriahzation are accompanied by shock waves in the baryon component. The duration 
of this process is close to the Hubble time tniz) |23|) so we restrict computations by tniz), 
which means the ending redshift Ze — 12 for the initial Zj = 20. One can expect that collisions 
of baryonic flows during the viriahzation of dark haloes result in considerable density variations. 
In order to study possible influence of such density variations on thermochemical evolution we 
conduct the calculations for a wide range of density in colliding flows: from the background to 
the virial value. As a characteristic value we adopt the virial density, 187r^ni,(l + z)^ (see e.g. 
j24j ) . where n?, is the background baryon density today, and consider the dependence of the 
thermochemical evolution on the initial gas density. Gas before the shock is assumed to be cold 
compared to the gas just after the front, so we consider strong shock waves. 

HD molecules form efficiently at low teperature in the presence of sufficient fraction of 
molecular hydrogen through the reaction 

D+ + H2 — >HD + H+. (6) 

For this reason we briefly discuss formation and destruction of H2 molecules behind shock fronts 

[ini mini 113. 

In the primordial gas II2 forms in interactions of neutral hydrogen with H" and Hg" ions effi- 
ciently born at high temperature. As shock waves significantly increase temperature it enhances 
formation of II2 and increases its abundunce ^j. It is known that HI cooling becomes inefficient 
at temperature ^ 10^ K, and in the primordial gas only the molecular hydrogen can provide 
further cooling to lower temperature. One can estimate a minimum fraction of H2 needed for 
effective cooling: the cooling time must be shorter then the Hubble time which condition fulfills 
when H2 fraction becomes greater than the critical value = 5 x 10~^ [211; increase of H2 
abundance shortens the cooling time. In a shocked gas at temperature ,^ 8 x 10'^ K collisions 
increase fractional ionization, and this further enhances formation of H~, H^ and H2. Figures 
1-3 show the evolutionary paths of the thermochemical state of a gas element behind the shock: 
the evolution begins at high temperature and follows a monotonous cooling accompanied by 
grownig H2 and HD fractions. It is seen that at high collisional velocities the electron fraction 
significantly increases at initial stages, which stimulates formation of H^ ions and H2 molecules 
and a strong decrease of temperature. The H2 abundance grows rapidly in the temperature 
range T ~ 10'^ — 10^ K, while at lower temperatures formation of H2 practically exhausts |25j . 
Already formed H2 molecules provide further cooling down to T ~ 200 — 400 K depending on the 
exact value of their abundance. If gas temperature falls down to T ^ 150 K deuterium rapidly 
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converts to HD molecules due to chemical fractionation (Fig. 4.). 

Fig. 5 presents the dependence of electron, H2 and HD fractions on the fluence rj = J ndt 
for several values of the shock velocity. Note that for the shock waves with Vg > S.Sa""*^ km s^"*^ 
the final abundance of H2 is greater than the limit 5 x 10~^, and as a cosequence the gas 
behind the shock can lose significant fraction of its thermal energy in one Hubble time. Since 
larger shock velocity corresponds to higher temperature, the maximum density behind the shock 
increases with the velocity p oc v'^ (eqs. 0]and|HI). In Fig. 5 this corresponds to larger fluence. 
It is obvious, that the characteristic time of the thermal evolution decreases as oc v~'^. For 
Vg — 4.6a^^ km s^^ H2 fraction equals ~ 7 x 10~^ and becomes sufficient for cooling down to 
T ~ 130 K; HD fraction at these conditions is 4 x 10^^. Further increase of the shock velocity 
results in an increase of molecular fractions and a decrease of temperature on a shorter time, 
with the dominant cooling provided by HD molecules (Fig. 6). Equating the cooling rates from 
H2 and HD one can estimate the critical temperature below which the contribution of HD into 
thermodynamics becomes dominant. For the adopted cooling functions it occurs at ~ 130 K. 
For the shock wave velocity ~ 4.6a^^ km s^^ gas cools down to this limit, and for larger velocities 
temperature falls below Tcr where the HD cooling dominates. 

In the velocity range 4.6a~^ ^ Vg ^ 8.7a~^ km s~^ the initial temperature behind the shock 
is insufficient for collisional ionization, and in the subsequent evolution the electron fraction can 
only decrease (Figs.l and 5). At such conditions H2 and HD abundances increase with velocity 
only because the reaction rates grow with temperature. As seen in Fig. 2 for shock velocities 
in this range the maximum fraction of H2 is approximately equal to x^^ ~ 8 x 10^'*. This is 
enough for an efficient formation of HD molecules and successful cooling down to T < 130 K. At 
V ^ 7a~^ km the HD fraction becomes sufficient for cooling down to the CMB temperature 
TcMB = 2.7(1 + z): due to the strong emission in rotational lines of HD gas temperature falls 
to several tens, ~ 30 K, approaching Tcmb(-2) at a given redshift (Fig.5). This is because 
HD molecules provide an efficient exchange of energy between the CMB and baryons through 
absorption of CMB photons and subsequent collisional de-exitation ^ El [^Hj • It is obvios that 
similar picture is valid at all redshift, and the final temperature of cold baryons is Tmin — ^cmb- 

For velocities Vs ~ 9.2a^^ km s~^ gas temperature behind the shock becomes greater than 
10^ K, which results in an increase of fractional ionization immediately after the gas ele- 
ment crosses the shock front: for = 10.4a~^ km s~^ Xg increases by factor of 2, and for 
Vg ~ 11.6a~^ km more than an order (Fig.l and 5). At such conditions chemical kinet- 
ics changes qualitatively - the enhancement of H2 formation is caused in this case by the two 
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factors: increasing reaction rates and a higher ionization fraction, resulting in more frequent 
catalystic processes H+e — >^H~^, H~^+H— >^H2 + e. The corresponding evolution of Xe{t), Xfi2{t) 
and XHD{t) looks qualitively different as seen in Fig.5: while for Vg < 9.2a''^ km s~^ an increase 
of velocity by 1 km s~^ produces insignificant changes in XH2(i) and XHD{t), for greater veloci- 
ties such an increase results in a considerable (half order of magnutude) increase of xb.2 {t) and 
xnD{t)- Thus, XHiit) and xhd(*) fractions at equal temepratures are higher for larger velocities 
(Figs. 2 and 3). For lower velocities the evolutionary paths xii2{T) and X}id{T) for different Vg 
practically coincide at temperature T < 10"^ K. One should stress, that in the considered range 
of velocities H2 molecules form primarely through H~ ions, the contribution from H2' ions is as 
a rule negligible. 

Thus, for velocities v > 4.6a~^ km HD molecules behind shock waves provide lower 
temperatures than H2 can do. It is worth noting that at z ~ 45 the CMB temperature is higher 
than the critical value Tcr at which HD cooling dominates. Therefore, under these conditions 
HD molecules can only heat gas, and at larger redshifts become unimportant. 

Everywhere above we adopted gas density in colliding flows equal to the virial value at the 
corresponding redshift. However, one can assume that in the process of merging of haloes a 
fraction of baryonic mass can be lost. In the intergalactic medium such "separated" baryonic 
flows can greatly expand, and their final density depends on collision velocity Vc, masses of the 
merging haloes, details of separation and so on. Subsequently such baryonic clumps can collide 
with gaseous components of other haloes or with each other. In these conditions the thermal 
evolution differs from that of denser baryonic flows. Let us consider how the thermochemical 
evolution depends on the density. Fig. 7 shows the temperature and the HD fraction versus the 
density for several shock velocities at two redshifts. It is clearly seen that at 2: = 20 and in the 
low velocity range v < 5. 80;^^ km s~^ only for densities close to the virial value gas temperature 
drops below the critical value where contribution from HD dominates. However, for higher 
velocities HD cooling remains efficient even for densities of one order of magnitude lower than 
the virial value. At v ~ 5.8a~^ km and the density close to the virial value only ~ 0.25 
of deuterium converts in HD, however it becomes sufficient to cool the gas down to the CMB 
temperature. For higher velocities this can occur for several times lower densities than the virial 
value. Collisions with higher velocities v > 11.6a~^ km s^^ change H2 kinetics: formation of 
H2 becomes more efficient due to significant increase of fractional ionization behind the front, 
and behaves similar to the collisions with virial density at Vg ^ 10.4a~^ km s~^ shown in Figs. 
1 and 5. These features are seen in Fig.7. In general, one can conclude that HD molecules can 
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also play a significant role in cooling of baryonic flows of low density. 



4 Discussion: formation of protostellar fragments 

Birth of stars is always accompanied with shock waves. This is unconditionally true for the 
first stars in the universe. The origin of shock waves can be connected both with the formation 
of the first protogalaxies in merging flows, and with supernovae explosions in already formed 
galaxies. Formation of the first protogalaxies implies separation of highly overcritical density 
perturbations from the Hubble expansion and a predominantly one-dimensional compression 
[261 127j . This is a source of shock waves in baryonic component. Similar processes take place 
in the hierarchical scenario of structure formation, where massive objects form in collisions and 
following mergings of less massive haloes. These processes can be treated as collisions of the gas 
and dark matter flows. A collisionless dark matter reaches the virial state apparently through 
the violent relaxation, while in gas component shock waves form. The shock wave velocity 
depends on mass of a forming protogalaxy: the velocity amplitude in a perturbation of mass M 
is close to the value 

Vc = \/3cr, (7) 
where a is the one-dimensional velocity dispersion |28j 

(8) 



, GM 

r = . 

2R 



Thus, the parameters of the shock waves are determined by the total mass of matter involved 
in motion, by the redshift at which the object forms, and so on. The efficiency of HD formation 
is sensitive to these parameters. Therefore, one can expect that the characteristcs of stellar 
population vary in galaxies of different mass. As the maximum abundance of HD depends on 
the shock velocity, and HD molecules cool the gas to much lower temperature than H2 molecules 
do, a typical mass of protostellar molecular clouds is expected to decrease with increasing mass 
of a forming galaxy |29j . 

Shock waves in the epoch of galaxy formation can be connected with explosions of the 
first supernovae. In these events much more powerful shock waves form: typical velocities 
can be greater ^ 100 km s~^, the corresponding temperature behind the front is ^ 2.8 x 
10^ K. At radiative stages when the gas temperature reaches ^ 10^ K, the H2 fraction becomes 
sufficiently high ^ 5 x 10~^ jl2l I30j due to high ionization fraction at the preceding stages. 
As a consequence, the gas temperature definitely falls to the lower values at which the cooling 
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is essentially determined by HD molecules. In these conditions fragmentation of an expanding 
shell can occur pT] . 

Due to isobaric compression the gas density behind the shock increases considerably com- 
pared to the initial value. For instance, for the velocity Vc > 7a^^ km s^^ the initial temperature 
is T > 5.8 X 10"^ K @, and when the gas cools down to Tcmb — 2.7(1 + z), its density increases 
more than 200 times. At such conditions fragmentation and formation of stars become pos- 
sible behind the shock j32j . Gravitationally unstable fragments can give rise to protostars or 
protostellar clusters. It follows that when cooling is determined by HD molecules, formation 
of low mass protostellar clouds becomes possible jHSlEl]- Indee, under these conditions the 
gas temperature falls down to ~ 2.7(1 + z) K which at z = 20 is 4 times smaller than can 
be provided only by H2 molecules. The Jeans mass Mj ~ SOT^^'^u^^^'^Mq behind the front is 
Mj ~ IST^riQ ^^"^Tq ^^"^Mq, where no, the gas density in a flow (a cloud) before collision, Tq is 
the temperature at the shock; here density at the front is taken 4no as for strong shock waves. 
If we assume that gas density in flows (clouds) before collision is equal to the virial value (see 
e.g. [24]), no = 187r^nfe(l + z)^, then for velocities Vc > 7a^^ km s^"*^ (or initial temperature 
To > 5.8 X 10^ K) the Jeans mass IS 

. . ..^M. (i±£)- . ... . :„3M. (A,)- (1±£)- 

where Sc is the ratio of gas density before collision to the background baryonic density. At the 
same time, when cooling is determined only by II2 molecules the typical gas temperature is of 
~ 200 K, and the corresponding Jeans mass ~ 13.5[(1 + z)/20]^^ times exceeds the value given 
by © 

- '-^ ^ iwhwT ^ '"'"^^ ir^r (1^)'" (^)'" 

In other words, the question of how massive are the fragments formed behind shock fronts 
depends on whether the cooling is determined by H2 or HD molecules. The density of the 
fragments is > 10 — 300 cm~^ depending on the redshift and the initial temperature. Subsequent 
collapse is isothermal until the fragment becomes opaque in H2 and HD lines, which takes place 
at the density ~ 10^ — 10^'^ cm At this stage, if the cooling is dominated by HD molecules, 
the Jeans mass is Mj ~ 3OTcmb^"^^^^^0 ~ 10^^(1 + zf^^MQ however when HD is 

underabundant the Jeans mass can be 2-3 orders greater |35j . Further evolution is determined 
by the accretion of gas onto the central core [361 137j . If the accretion rate is below the Eddigton 
limit, the mass of a forming star is comparable to the initial mass of a protostellar cloud, in the 
opposite case it can be much lower |351 136j . Thus, one can expect that a typical mass of the 
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first stars born in protogalaxies of liigher masses (corresponding to liigher collisional velocities) 
is shifted towards the lower end due to cooling by HD molecules. 

It is readily seen that since the overall thermal evolution of low density flows differs from 
that of higher densities, the final value of the Jeans mass and its dependence on redshift will 
differ from the above value. Fig. 7 shows the Jeans mass versus the gas density in the flow. In 
the range of collisional velocities 5.8a~^ < "Uc < 8.6a^^ km s^^ only flows with the initial density 
very close to the virial can have Jeans mass Mj smaller than W^Mq, which may correspond 
to the mass of a protostellar cloud. However for v > 8.6a~^ km s~^ the Jeans mass becomes 
< IO'^Mq for the initial density 4 times lower than the virial. As mentioned above, for higher 
velocities gas cools down to the CMB temperature, and the Jeans mass is < 10^ Mq, what is 
seen also from Q - on Fig. 7 flat parts of the lines for the velocity > 8.6a^^ km s~^ reflects 
this circumstance. Thus, the Jeans mass is considerably higher than Q only for flows with a 
low collisional velocity and a low density. High-velocity collisions Vc > 8.6a^^ km s^^ provide 
cooling down to the temperature T ~ Tcmb even for low density flows. 

Let us consider collision of flows whose density is equal to the background value, i.e. (5c — 1, 
p/Pvir ~ 6x10^3. It is seen from Fig. 7 that even for high-velocity collisions Vc = 11. 6q ^ km s ^ 
gas cannot cool sufficiently in one Hubble time, and the Jeans mass is quite high: ~ 10*^ -W Mq. 
Moreover, the free-fall time for such low densities is greater than the comoving Hubble time. 
Under these conditions baryonic objects cannot be formed. However, further increase of the 
collisional velocity, Vc > 11. 6a^^ km s~^, makes the gas behind the shock able to cool down to 
the temperature < 1000 K during one Hubble time. For instance, for low density flows collided 
with the velocity Vc — 19. 2q^^ km s^^ the final temperature behind the shock is ~ 200 K, and 
can reach lower values for higher velocities. Under these conditions HD molecules will dominate 
in gas colling, and the Jeans mass becomes as small as ~ 7 x 10^ Mq. 

5 Conclusions 

The infiuence of HD molecules on the thermochemical evolution of the primordial gas behind 
shock waves possibly formed during the epoch of galaxy formation has been studied. 

1. We showed that deuterium converts efficiently to HD molecules and the contribution of 
HD to coolig becomes dominant for the shock waves with velocities ^ 4.6a~^ km s~^ 
(a ~ 0.5). Behind such shock waves the conditions are favourable for fragmentation and, 
as a consequence, for formation of protostellar clusters. 
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2. For shock velocities ^ 7a km s gas is able to cool down to the CMB temperature. 
Under these conditions Jeans mass depends only on the redshift and the initial density: 
Mj < 2.4 X 10^Mq{1 + z)°•^(To(5c)~°•^ for virial haloes {Sc = IStt^) at z = 20 this corre- 
sponds to Mj < 10^ Mq. 

3. At z ^ 45 the CMB temperature is close to the critical value T^^, at which the contribution 
from HD molecules to the total cooling is comparable to that from H2 . At these conditions 
HD molecules begin to heat gas, and at higher redshifts become unimportant in thermal 
history of baryons. 

4. For densities of colliding flows smaller than the virial value the efficiency of HD molecule 
formation decreases. In particular, at 2; = 20 gas temperature behind the shock with 
V ~ b.Sa"^ km s~^ drops substantially only for a density close to the virial value. However, 
for the shock velocities ~ 8.6q:~^ km s~^ HD molecules are important in cooling for 
densities of 2-3 times lower than the virial value. For v ~ 11.6q;~^ km s~^ HD cooling is 
effective for densities close to the background value, and almost all deuterium converts to 
HD for collisions of less dense than if the flows were virial. For the gas density equal to the 
background value and for the velocity v ^ 19.2a~^ km s~^ temperature drops to < 200 K 
and HD molecules begin to dominate radiative cooling. 
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the temperature range of interest (T < 10^ K) this function coincides with that given in [18], 
while at higher temperatures, where Lipovka et al predict an order of magnitude enhanced HD 
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Figure 1: The evolutionary paths connecting fractional ionization and gas temperature behind 
shocks in collisions with velocities 2.3, 3.5, 4.6, 5.8, 7, 8, 9.2, 10.4, 11.6 km - from top 
to bottom; initial points of the evolution are marked by filled circles. 
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Figure 2: The evolutionary paths connecting variations of temperature and H2 concentration 
for the same velocities as in Fig. 1 (left to right). The vertial line depicts the critical value 
5 X lO""* needed for baryons to cool in one comoving Hubble time |24j . 



14 



Ti| — I I I — I I I — I I I — I I I — r 




10"^ 10"^ 10"^ 10'^ 10"^ 



x(HD) 

Figure 3: The evolutionary paths for the abundance of HD molecules (see Fig. 1). 
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Figure 4: The relative contribution of radiative losses in the lines of H2 and HD into the total 
cooling (see Fig.l). Squares mark the initial points, circles - the final points of the evoluton for 
the velocities 2.3, 3.5, 4.6, 5.8, 7, 8, 9.2, 10.4, 11.6 a~'^ km s^^ 
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Figure 5: Ionization, H2, HD fractions and temperature behind the shock front versus 77 for the 
velocities 2.3, 3.5, 4.6, 5.8, 7, 8, 9.2, 10.4, 11.6 a'^ km s""*^. Circles on temperature lines mark 
the final points of the evolution. 
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Figure 6: The relative contribution of radiative losses in the lines of H2 and HD into the total 
cooling (see Fig. 5). Filled (H2) and open (HD) circles mark final points of the evoluton for 
velocities 2.3, 3.5, 4.6, 5.8, 7, 8, 9.2, 10.4, 11.6 a"^ km s'^ (left to right). 
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Figure 7: Temperature, HD fraction and Jeans mass reached behind the shock fronts with the 
velocities 2.9, 5.8, 8.6, 11.6 km s"^ for two redshifts: z = 20 - solid, z = 15 - dashed lines. 
Shaded area on the top panel shows the region where the HD cooling prevails. The density in 
units of the virial value for z = 20 is shown at the upper axis. 
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